SEMIANNUAL REPORT FOR 


/P_^_ tf)p< 

Z/4 c-l'-roo 


Grant No. NAG 1-500 



SEMICONDUCTOR SUPERLATTICE PHOTODETECTORS 


July 1, 1985 to December 31. 1985 


Submitted to 


Dr . Ivan Clark 


National Aeronautics and Space Administration 
Langley Research Center, Mail Stop 473 
Hampton, VA 23665 


<kiSi-CB> 180300) SEHICCNDUCTC 
f EOIODE1ECIOBS Setiannaal Sepc 
31 Dec. 1985 (Illinois l/riv.) 
K3IS HC A03/HJP kQ 1 ■ 


S0PEHLA1TICE 
, 1 JUX. - 
*1 P Avail: 
CSCI 20L 


G3/76 


N87-27520 


Unclas 

0030524 


Prepared by 


S. L. Chuang 

Department of Electrical and Computer Engineering 
University of Illinois 
Urbana, IL 61801 

K. Hess. J. J. Coleman and J. P. Leburton 
Department of Electrical and Computer Engineering 
and Coordinated Science Laboratory 
University of Illinois 
Urbana, IL 61801 


SEMIANNUAL REPORT FOR 


Grant No. NAG 1-500 


SEMICONDUCTOR SUPERLATTICE PHOTODETECTORS 


July 1, 1985 to December 31, 1985 


Submitted to 


Dr. Ivan Clark 


National Aeronautics and Space Administration 
Langley Research Center, Mail Stop 473 
Hampton, VA 23665 


Prepared by 


S. L. Chuang 

Department of Electrical and Computer Engineering 
University of Illinois 
Urbana, IL 61801 

K. Hess, J. J. Coleman and J. P. Leburton 
Department of Electrical and Computer Engineering 
and Coordinated Science Laboratory 
University of Illinois 
Urbana, IL 61801 


TABLE OF CONTENTS 

page 

I. INTRODUCTION 3 

1 . Period .......... 3 

2. Reporting Date . ....... 3 

3. Technical Personal ....... 3 

II. TECHNICAL PROGRESS 4 

III. FUTURE WORK 6 

IV. PUBLICATIONS 6 

V. TRAVEL 7 

APPENDIX: 


"Impact Ionization Across the Conduction Band Edge Discontinuity 
of Quantum-Well Heterostructures'' by S. L. Chuang and K. Hess.... 8 


3 


I. INTRODUCTION 

The research grant NAG 1-500 entitled "Semiconductor Superlattice 
Photodetectors" was awarded to the University of Illinois at Urbana-Champaign 
by National Aeronautics and Space Administration-Langley Research Center 
on June 27, 1984. The grant was continued on July 1, 1985. Dr. Ivan Clark 
is the Technical Officer, and Mr. John F. Royall is the Grants Officer. 

This report is the third semiannual report. 

1. Period: 

July 1, 1985 to December 31, 1985 

2. Reporting Date: 

January 10, 1986. 

3. Technical Personnel: 

S. L. Chuang Assistant Professor of Electrical and Computer 
Engineering 

J. J. Coleman Professor of Electrical and Computer Engineering 

K. Hess Professor of Electrical and Computer Engineering and 

Research Professor of Coordinated Science Laboratory 
J. P. Leburton Assistant Professor of Electrical and Computer 
Engineering and Research Assistant Professor of 
Coordinated Science Laboratory 


Two research assistants 


II. TECHNICAL PROGRESS 


During the past half-year period, we continued the work on the new 
superlattice photomultiplier with the following effects investigated. 

The detailed analysis is included in the Appendix of this report. Here 
we give a brief description. 

(1) The effect of multiple subbands 

Since the quantum-well regions may contain several subbands, each of 
which may be occupied by electrons depending on the doping concentrations, 
it is important to include the multi-subbands in calculating the impact 
ionization rate. The electrons occupying the higher subbands require a 
smaller amount of energy to get out of the quantum well; thus, those higher 
level subband electrons contribute significantly to the impact ionization 
rate. The results of the multiple subbands have been calculated. (Please 
see page 17 of the Appendix.) 

(2) The nonparabolicity effect of the band structure 

Since the electron energies involved are of the order of 0.3 eV or 
higher, the nonparabolic band structure may be important. That is, the 
effective mass of the electrons may be larger than the effective mass of 
the electrons at the band edge of the T valley. Thus, an averaging procedure 
has been used to find an average effective mass of the electrons that 
takes into account the nonparabolicity effect. (Page 26 of the Appendix) 


(3) The effect of the quantum-well size 

The quantum-well size determines the energy levels and the Fermi level 
for a given doping concentration. Changing the well width will change the 
impact ionization rate. It is found that for the same doping concentration 
in the quantum wells, the well with a thicker dimension will have a larger 
ionization rate because there are more electrons available for ionization 
in a thicker well . 

(4) The effect of the band-edge discontinuity 

The band-edge discontinuity affects the threshold energy for impact 
ionization. A smaller band edge discontinuity will give a larger impact 
ionization rate because the threshold energy required for impact ionization 
is smaller. (Page 18 of the Appendix) 

(5) The effect of doping 

The doping concentration decides the number of electrons available 
for impact ionization. Increasing the doping concentrations will increase 
the impact ionization rate. The results of various dopings are shown in 
Figure 7 of the Appendix and discussed on page 19 of the Appendix. 


III. FUTURE WORK 


Recently, Capasso at AT&T Bell Laboratories has independently 
performed experiments on similar types of structures and obtained 
encouraging evidence for impact ionization of the kind discussed in 
this research. These experiments may call for an extension of our theory 
to include effects such as the tunneling-supported impact ionization and 
the impact ionization from deep impurity centers. We shall investigate 
these effects and others such as the thermionic emissions. These studies 
will certainly provide more insight into the physics of the superlattice 
photodetectors . 


IV. PUBLICATIONS 

The following manuscripts submitted for publication were supported 
either fully or partially by this Grant. The support by this Grant has 
been acknowledged in these manuscripts. 

1. S. L. Chuang and K. Hess, "Impact ionization across the conduction 
band edge discontinuity of quantum-well heterostructures", J. Appl . Phys . , 
accepted for publication. 

2. S. L. Chuang, "Lateral shift of an optical beam due to leaky surface 
plasmon excitations", J. Opt. Soc. Am., accepted for publication. 

3. S. L. Chuang, "Anisotropic wave propagation in a superlattice waveguide, 
submitted for publication. 


V. TRAVEL 


S. L. Chuang attended the 1985 International Electron Devices Meeting 
at Washington D. C., Dec. 1-4, 1985. He also attended a short course on 
"Digital III-V Device and Circuit Technology" during that meeting. 
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"Impact Ionization Across the Conduction Band Edge Discontinuity 
of Quantum-Well Heterostructures" By S. L. Chuang and K. Hess. 


IMPACT IONIZATION ACROSS THE CONDUCTION BAND EDGE 
DISCONTINUITY OF QUANTUM-WELL HETEROSTRUCTURES 


S. L. Chuang 

Department of Electrical and Computer Engineering 
University of Illinois at Urbana-Champaign , Urbana, IL 61801 


K. Hess 

Department of Electrical and Computer Engineering 
and Coordinated Science Laboratory 
University of Illinois at Urbana-Champaign, Urbana, IL 61801 


ABSTRACT 

Impact ionization across the band edge discontinuity of quantum-well 
heterostructures is studied theoretically. We consider a heterolayer structure 
of alternating Ai^Ga^^As and GaAs layers where the GaAs layers are heavily 
doped with donors. Thus a large number of electrons is confined to the quantum- 
well region. Incident electrons are heated up by applied electric fields and 
collide with the electrons confined in the well regions. Both the ionization 
rate as a function of the incident energy, and average ionization rates are 
computed. Device applications of such multiple quantum well structures and 
the possibility of a complete analog to the conventional photomultiplier are 


discussed 
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1 . INTRODUCTION 

Novel types of photodetectors using superlattice structures have been 
proposed recently 1 10 . These include the enhancement of impact ionization in 
a superlattice for avalanche photodiodes, 1 6 an infrared detector using free 

7 , 

carrier absorption in a heavily doped quantum-well structure, a photodetector 

O 

using the transient displacement current effect in a sawtooth superlattice, 

9 10 

strained-layer superlattice photodetectors ’ and doping (n-i-p-i) photodetec- 
11 12 

tors. * The operational principles of all these new devices are very inter- 
esting and additional studies are needed to characterize these photodetectors 
completely. 

In this paper, we consider a heterolayer structure which consists of 
alternating A£. x Ga^_ x As and GaAs layers where the GaAs layers are heavily doped 
with donors (Figure 1). The donor are electrons confined mostly to the quantum 
wells. Incident electrons (which for example can be generated optically) are 
accelerated by the external field and collide with the electrons confined in the 
wells as shown in Figure 1. This results in impact ionization of some of the 
electrons bound in the quantum wells. The wave functions for the incident 
electrons (three-dimensional) and the confined ‘electrons (quasi-two dimensional) 
are used to calculate the electron-electron interaction due to the Coulomb force 
with screening effects taken into account. The quantum mechanical transition 
rate for a system involving an initial state of one free electron above the well 
and one bound electron in the well and a final state of both electrons free is 
calculated using the Golden rule. The process is schematically displayed in 
Figure 2. The electron-electron interaction in the momentum space is illustrated 
in Figure 3 and will be explained in the next section. The average ionization 
rate over the momentum (energy) distribution of all the incident hot electrons 
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is then calculated. Potential device applications are discussed in a final 
section. 

Our formulation for the impact ionization for the quantum-well structure is 

13-17 

similar to the inverse Auger effect in bulk semiconductors or in a quantum- 

1 8 * 1 9 

well heterostructure. The new feature is that we are interested in impact 

ionization across the conduction band edge discontinuity while previous studies 
are considering transitions across the energy band gap. Thus, the threshold 
energy for impact ionization is smaller (of the order of conduction band edge 
discontinuity) than that of the interband process. The required bias voltage 
may be smaller since a smaller electric field than that of the conventional 
avalanche photodiode is needed. 
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2. ANALYTICAL CONSIDERATIONS 

In this section, we present the detailed formulation of the impact ioniza- 
tion of electrons across the conduction band edge discontinuity in a heterolayer 
structure. First the transition rate due to the electron-electron interaction 
is derived, then we discuss the energy distributions of the electrons confined 
in the quantum-well region and the incident hot electrons. Third we present the 
impact ionization rate and its numerical evaluation. 


2.1. Transition Rate 

Consider an incident free electron interacting with a bound electron in a 


quantum well; the total average transition rate per unit volume P Is given 


by 


19,20 


P tc -7 I I l l p! 1 ! 2 «*!> « k V t 1 - 

,+ y k, k„ 


( 1 ) 


h z S h 12 


k l k 2 

where f is the Fermi distribution and P is the quantum mechanical transition 
rate (per unit time) between the initial state of electron 1 in state k, , 

X 

electron 2 in state and the final state of electron 1 in state k^ and 
electron 2 in state k^ (Figure 2) due to the electron-electron interaction: 


k ^ t f 

P 1 2 = I 1 ! < 1 2 |H |l’2’>| 2 5(E + E - E - E) 

it i t “ S l Z l Z. 

12 


( 2 ) 


The square of the matrix element | < 1 2 | H |l'2'>| 2 is gi ven by 


o * + + e exp(-q|r - r |) * > > ,2 

I < 1 2 | H 1 1 * 2 ’ > | 2 = 4 | J (r ) (r ) — “ + v (r i ) MV dr i dr 2 1 

k l k 2 4 tt e ! r i - r 2 1 ^ ^ 


( 3 ) 
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for the interaction by a screened Coulomb potential. Here q is the screening 

parameter which is discussed in Appendix A. We have neglected the difference 

in the dielectric constants of the two materials which is a good approximation 

for A&GaAs-GaAs . The factor of 4 in Equation (3) takes account of effects of 

14 15 

the exchange and preservation of the spin ’ , considering the integrations 

over icj and ^ » an< * 


<V 2 |12> h K ( V V ( V 


*1 2 
exp(i^ ^r^) 

N 


’ /aE~ exp(i ^2t* ? 2t ) Sin I^ z 2 
v zn zn 


(4) 


<f 1 f 2 tl'2’ > = +,(?!> * ,(t 2 > 

k l 2 

exp(i^ 1 «r 1 > exp(iic 2 .r 2 ) 

SB • - - 

/V /V 


(5) 


where k^ , k^ , k 2 , r^ , r 2 are the general three-dimensional vectors and k 2j _, 
are two-dimensional vectors in the quantum well interface planes. The Bloch 
functions associated with Equations (4) and (5) have been ignored since the over- 
lap integral is assumed to be one for states in the conduction band. The effec- 
tive quantum well width L is defined such that 


Qn 


. mr 2 

2-* (L Z n 


( 6 ) 


where E^ n is the energy level of the nth subband, which is calculated numerically 

taking into account the finite barrier height AE , i.e., the conduction band 

c 


edge discontinuity (See Appendix B) . For convenience, we further define: 


(7) 
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. n tt 

4 + 4 A AAA n. 

k = k ± z — = xk +yk ±z 

2n 2t L x y L 

z zn 


r 2 " r 2t + zz 2 " XX 2 + yy 2 + ZZ 2 


(8) 


and replace the summations £ by integrations. One then obtains 

£ 


|< 12 | H s | l ' 2'>] 2 . (^) 2 (^) [|£i .‘^2 + q 2 )2 [ V k 1 - k V f 2 

+ 6 » .] 
k l k 2 ,k l +k 2 

and 


(9) 


4 , L 2 dk 
- — iLIL <- zn 


dit_ 


tr 2 
e 


rO OO) 


zn 


Here we have used the fact that f(jc^) and f(^) « 1» also f(ic^ t ) = f(ic^) = f(ic^) 
(which will be discussed later) and 


* ++ «(E, + E„ - E - E ) 

^ 1’ 2n ; " J .3 r ,* .2 2.2 

, , (2ir) [ |k - k j + q J 

k, >Ku 

lz— th 

* i 

k„ >k 
2z— th 


(ID 


+ 1 > >+ ■>* 
k =k +k~ -k, 
2 1 2n 1 


The lower limit in the integration for k is due to the fact that both k, and 

1 1 z 

* i 

k need to be larger than a certain threshold k , in order to get both electrons 
2z th 

f f 

1' and 2' out of the quantum well. The condition k. > k , imposes another con- 

2z — th 


V « « V V 

dition on k, because of the relation k = k ± k > k , . One finally 

lz 2z lz L lz — th 

zn 


obtains 
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, / . , , , , nit , 

k < k < k ± k , 

th — lz — lz L th 

zn 


( 12 ) 


One notes that the above condition for k, also holds for k„ obviously. 

lz 2 z 


An upper bound of the integral of Equation (11) without the constraints on 


kj z and k^ z has been obtained analytically for a parabolic band structure 


21 


s (it ,£* ) 

u 1 2 n 


dk. 


5(E 1 + E 2 - Ej - E 2 ) 


all ( 2 it ) 3 [|^ - £[ j 2 + q 2 ] 2 


m* 


1*1 - k 1ni 


8. 2 « 2 q 2 |Sj - it| n P + , 


.♦± .2 


£=£+£+ -it 
2 1 2 n 1 


(13) 


.+ 


The geometry pertinent to the evaluation of S(k^ ,k 2 ) and its upper bound 

S (&, ,£_) is shown in Figure 3, for + z 7 ^— (n=l). Due to energy 

u 1 2 Z zn 2 t L 

t zn , 

^ .y 1 

and momentum-conservations, the four vectors , k^ * k^ and k^ lie on a spherical 

-► 2 2 

surface assuming that the E-k relations are parabolic, for all four vectors 

» » + + 

Due to the above mentioned constraints, on both k, and k , the S(k ,k ) 

lz Zz 12 

integral will be over the unshaded spherical surface in Figure 3 for which 

* t _ » + + 

k .< k, < k, + k , , while the integral for S (k ,k ) is over the whole 

th— lz— lzL, th ul2 

zl 

spherical surface. The latter integration can be peformed analytically by a 
coordinate transformation to the polar axis ^ ^ , and results in Equation (13) 

2.2. Electron Distribution Function f(ie 2 ) t * ie Quantum-Well Region 


+±. 


The electron distribution function in the quantum-well region f(k 2 ) or f(k 2 ^) 
is taken as 

r (n) 


f ( J 2t ) - ( 


1 if 0 < k 2 t - K 2 t 
0 otherwise 


(14) 


where 
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JP> 

zt 


/, * 2 

-2 E F 

■n zn 


( 15 ) 


is the maximum of k^ t for electrons in the n-th subband. 
The electron concentration is given by 


n = / f (E) g(E) dE 
E 

Ql 


= I Hr - V ln[l + e ( V E Qn )/k B T ] 

n irfi L B 
z 


(16) 


where E is the n-th quantized energy level of electron 2 as given in Equation 

(6). Defining a characteristic concentration n , 

c 


_ • V 

nC Trfi 2 L 


(17) 


We have for E„ - E„ » 

F Qn B 1 ’ 


(E_ - E )/k T . 
n u , F Qn B 

c occupied 


— = l 

OC 

subbands n 

From n » N n and the above relation or (16), we obtain Ej” anc * t ^ lus ^ 2t* 


(18) 


2.3. Electron Distribution Function f(k^) for the Incident Hot Electrons 


We assume that the incident hot electrons in the Ail Ga As region (region a) 

X 1 ~“X 

is Maxwellian in shape, i*e. , the spherically symmetrical part is given by the 
expression: 

_j*2 3/2 _R )/k T 

(19) 


. , , irf . 2 3/2 -E^CSj/kT 

f(k ) = 4n (— 5 ) e 0 0 B e 

v y o i m 


2m k n T 
a B e 

* 

with an electron temperature T , where ra^ is the effective mass of the electron 

in the conduction band in the A£ Ga As region and the index 0 characterizes 

x 1 -x 
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the electrons in this region. Since the energy 

.2 


E 0 = 




2m 


* (ki. + kL + kj_) 


Oy Oz' 


( 20 ) 


is measured from the conduction band edge in the A l Ga As region, and 

X 1 —*x •» 


k. = k. x + k y + k. z , 
0 Ox Oy Oz 


( 21 ) 


the energy E^ and k^ vector in the GaAs region (region b) are related to and 
by the equations 


k, = k 
lx Ox 


( 22 ) 


k = k 
ly Oy 


'5-* kj: + 4E - 2-_ k 2 

2m 0 c 2m* 1 

a b 


(23) 

(24) 


where m^ is the average effective mass of electrons in the T-valley of the GaAs 

•k 

region. (We use the average effective mass m^ to include the nonparabolicity 

* 

effect. The effective mass m used in the previous equations is replaced by 
m fe in the computations. See Appendix C) . Instead of integrating over dk^ , we 
integrate over dic^, 

dl 


dk dk 

/ 3 f(k,) * / J ^n^O^ 

(2ir) J 1 (2X)- 3 0 ° 


2.4. Average Impact Ionization Rate 

The average ionization rate <■— > is defined as the total transition rate per 
incident electron and therefore: 


10 - 


<-> 

T 


tr 


n 


0 


= / 


- / 


dk„ 


irfi 


2 3/2 


2 2 2 


’ exp[ '^T <k o* + V k o 2 >l 


( 2u) 3 

H * ) 

2m k_T 
a B e 

* 

2m k T 
a 6 i 

dk 0z , 

J2 1/2 

7m ^ 

2 2 
-n k 

„ r Oz 

2ir 1 * > 

^ 2m k T 

a.B e 

ex PL~ * 

2m k T 
a B e 


T<E Oz> 


(25) 


where 


4 L 2 
e 4 if /• zn 


dl 


x(E n ) 
Oz 


o (-f 1 ) / «*2t» s<k V k ~„> + 3(1?, A- )| 

e h (2ir) L 


zn 


(26) 


which is the ionization rate for an incident hot electron with the energy 

E ° ‘ E °* ' 17 k °* 

a 

and 


(27) 


k n = k o = 0 
Ox Oy 


(28) 


Since the integrand in (25) is sharply peaked near k = 


Ox 


the approximations 

^7 


lCQx^O T ^ E 0z^ 

k0y=O 


k. ~ 0, we have used 
Oy 


(29) 


and carried out the dk^ and dk^ integrations in (25). By doing so, the inci- 
dent electrons are essentially z-directed. 


* 

k- " z k_ 

0 Oz 


( 30 ) 


Thus, for the transmitted electron we have 


11 - 


k - z k 
1 lz 


(31) 


and 


r k!T + AE 


2m 


Oz 




2m, 


lz 


(32) 


The average ionization rate <— > is obtained by integrating Eq. (25) using Eq. 

(26). An upper bound of <-^> , <~> u » is obtained using the upper bound of ^ ^ ♦ 

Dropping the subscript z and using E rt instead of E. , one obtains for the inte- 

0 Oz 

grand of Eq. (25) 

, 4 , L die 

tcOu 3 % ff ^ / 77 jr- ^ 1 V k V%„> ♦ • <»> 

0 e "n L 

zn 

Below we first discuss the exact numerical integrations of -y and <-> in Eq. 

1 ® 

(26) and Eq. (25). Then formulas for upper bounds, ">VT '» which lead to analy- 

x T(E o ;u 

tical expressions for <— >^ are presented. 
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3. THE NUMERICAL INTEGRATIONS 

3.1. Integrations and <— > from Eq. (26) and Eq. (25) 

To obtain -A - ~r , we need to evaluate the electron-electron interaction 

I± 0> 

integrals S(k^ ^k^) F° r tbe following discussion, we use 

k = let • The integral S(& ,^T) is obtained similarly, replacing let by ic7 , 
^ 1 / 2n 2n 

° r by TT^~ # For a P ara bolic band structure, and using 
zn zn 






it + 

zk 

it 

lz 

it' + 

* t 

zk , 

2t 

2z 

+ E„ - 

? 

■ E - 

2 

1 


-j 6[(k’ - k ,0 )(k' - k' B )] 

^2 lz 1 z lz lz 


where the two roots for k, , are 

1 z 


k. 

, 1 z L 

k a = 5£ + 

lz ' 2 


mr 

L 2 
zn ^ 


) + it' .it - k' 2 

; it 2t it 


k + - 2 - 
'6 _ 12 L zn 


lz L 2 
zn ^ 


) + k lt- k 2t - 


Therefore, the k integration can be performed analytically using the properties 
i z 

t 

of the delta functions, and k is restricted as given in Eq. (12). The other 

J. z 

integrations can be simplified as 
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(n) 


I d *2t I <t ' L “ k 2t dk 2t C C “it dk it 2 ’ ! d *u 


2 IT 


(39) 


where the two angular variables $ 2t for and <|> lt for k have been trans- 
formed to only one angular variable <|>^ 2 which is the angle between k^ t and k^, 
with a 2it factor added as a result of this transformation because the integrand 

I V 

depends on <^ 2 only. The upper limit of k^ variable, K^, f s obtained from 

* a 1 S 

the condition that k, and k, in Eq. (37) and Eq. (38) should be real, i.e., 

1 z 1 z 


lz 


nir 

L 2 
zn 


) + * 


, k„ cos <fc,„ 
It 2t y 12 


k iU° 


(40) 


This condition is easily included in the numerical code* If it is violated, 
f a f S 

k and k become imaginary and the integrand vanishes because the variable of 
lz lz 

» 

integration, k. , is real. Explicit solutions of Eq. (40) can also be used, 

1 z 

f 

which give the upper limit for k^. 


-k / k - — 

’+ _ 2t A / lz L 2 k. 2 

K it - — + J ( — + (-ft) . 


(41) 


Then 


*<V 


can be written as: 


i e 4 m* L 2 , 2ir 

w 1 ■ -rf ^ ** dk - d ^ 

zn 

k’ + 2 it 

U 0 U k it dk it f 0 d *12 F+ < k U’ k 2f>12> 


» _ 


__ I 


+ k U dk it l ' d +12 F ' (k lt' k 2fh2» 


(42) 
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where 


F (k lt’ k 2t ,4> 12 ) = 


1 


nu 


— ** 11 "1 9 

of k, - L • 

2 { lz zn + k. 

L 2 J 1C 


k 2t COS +19 " k ^ ^ /2 


{ — 


1 


12 It 
1 


2 L . ' cu 2 , 2.2 ri ’2 „ . '0.2 2,2 

[ k l t + (k lz " k lz > + q ] [k 1(; + (k lz k lz ) + q ] 


} 


(43) 


f <x T 8 T nir f 

for and between k^ an d k^ + and zero for the appropriate 

zn 

_ » i _ 

terms otherwise. F (k, ,k , <(>, „) and K lt are obtained by replacing mr/L by 

it Zt 1Z It zn 

+ ? + 

-nTr/L in the F function as well as in R. from Equations (41) and (43). 
zn It 


The last integration f dk in <— > must be performed from a minimum k 


Oz 

which is obtained from Eq. (12) 


Oz ,min 


k, > 2k . + 7-S- 
lz — th L 

zn 


(44) 


and the energy conservation relation Eq. (32), 


* 

2 2m 


Oz ,min 


-* 1 ?-) - 
m^ zn ti 


AE 

2 c 


(45) 


3.2. The Upper Bounds of the Integrals and <~> u using Eq. (33) and Eq. (25) 


The analytical result for the upper bound of 


t(E 0 ) 


is 


4 * _ 0 

e m, L 2 

_ b / znx 

x(E n )u ~ 2J U J 


1 


0' 


e ff 


4ir^qL 

zn 


[G + G ] = 


T< V 


u t(E 0 ) 


(46) 


-IS- 


C' = 


/ 

/ 


K 2t 2 

Hr) + (■ 


k + 2 *- 
lz L 2 


zn 


k , + 

1 2 L 


) -(-7 — -) 


-tan’ 1 / (^) 2 + ( 

v q ' ^ q 


k + 
lz L 2 
zn 


k, *21 
, lz L 

) + tan_1 ( — 


Here the G term is due to the term with k„ , and the G„ term due to that 

2n 2 

£_ in Eq. (33), and 
2n 


i 00 

<i> - 1. * 


dk 


Hz. 


2ir 


Oz ,min 


dk 

°° Oz 

+ ; k- ^ 

Oz ,min 


, irfi 2 ,1/2 , * 2 ko z 

l — 5T “) eit P ( * 

W. 2m a k B T a 

-2 1/2 H 2 k? 

(-#— ) a*P (- ° 2 

2m k T 
a B e 


2ra*k T 
a B e 


t( e o> 


t(, W 


u 


(47) 


with 


(48) 


which can be integrated numerically by Simpson's rule 
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tCEq) 


, and the 


4. NUMERICAL RESULTS AND DISCUSSIONS 

The direct numerical integrations of the ionization rate 
average ionization rate <— > in Eq. (25) is straightforward. The upper bounds of 

T 

the ionization rate and the exact numerically integrated results as functions of 
the energy of the incident hot electrons are shown in Figures 4(a) and 4(b). One 
can see that the ionization rate increases rapidly when the energy E^ is above 
the threshold value E^ 


E i = 


+2 


* 2 

^( 2 *, 


nir 


„ * Oz.min . *'• th L 
2m ’ 2m, zn 

a b 


) - AE 

' r 


(49) 


for each subband n * 1, 2, 3. The overall minimum threshold energy for impact 

ionization will be due to the highest occupied subband of course. For a quantum 

well size L ** 200A, A l mole fraction x * 0.35 in the A % Ga, As region, and N = 
z x 1-x D 

18 -3 

5 x 10 cm , there are three occupied subbands among the six quantized levels 
in the well. The Fermi level is 0.1092 eV measured from the conduction band 

t I 

edge is the GaAs region. The threshold value for k , (the minimum value of k 

t n 1 z 

or k„ for electrons 1' and 2' to get out of the quantum well) is given by 
2z 



/ 


/2m *( AE - 6), 


(50) 


where 6 accounts for the tilting of the quantum well due to the large applied 
field (Fig. 1). Without space charge effects we have 


6 = eFL (51) 

z 

where e is the magnitude of the electron charge. In all our computations we 
assume a field strength F = 10 kV/cm and a corresponding electron temperature 
of 2300° K at a lattice temperature of T * 77° K. (No accurate data for in 
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A £. Ga As are available to our knowledge. The above value is a result of Monte- 

X x —x 

Carlo simulations in GaAs. In reference 23 there are some discussions on F-T 

e 

relations in AJ^Ga^^As.) The necessary parameters used here are taken from 
references 24 and 25. 

The average ionization rate <X> involves the integration over the energy 

Eq and the distribution function. The upper bound of the impact ionization rate 

x(E " )u as a ^ unct ^ on t ' ie energy Eq is shown in Figure 4(a) for each subband 

and also for all subbands combined. These curves are much easier and faster to 

generate than the results of the exact numerical integration which are shown in 

Figure 4(b) since the results in Eq. (46) are analytical expressions. They also 

can serve as a means of checking the numerical accuracy of the exact numerical 
* 

integral in Eq. (42). This is because ~ qT y u from Eq. (46) is the upper bound 

1 i 

of the numerical integral of — ^ ^ in Eq. (42). is the integration over 

0 0 +> + » 

over the whole spherical surface in Figure 3 , for which all k^ and k^ of the two 
electrons are allowed everywhere- on the spherical surface. ^ is the inte- 

gration over only a portion of the spherical surface for which the condition 

k\ < k, < k, + k \ is satisfied. If one relaxes this condition in 

th — lz — lz L th 

zn 

Eq. (43), the results of integration (42) will be over the whole spherical 
surface, which is the upper bound. Thus a numerical test in Eq. (42) by simply 

I 

removing the constraints on k. can be checked exactly with the analytical 

lz 

expression in Eq. (46). This process determines the number of discretized 
intervals for each integration variable in Eq. (42). One then adds the con- 

i i 

straints on both k, and k„ to obtain the exact numerical integrals. 

lz 2z s 


4.1 . The Effect of Multiple Subbands 

As shown in Figures 4(a) and 4(b), the ionization rates due to higher level 
subbands are clearly very important because electrons in those subbands are in 
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higher energy levels and thus easier ionized. Figure 4(a) shows the upper bounds 

of the ionization rates for each subband and the total contribution is shown as 

the solid curve. The total rate in Figure 4(a) has a discontinuity whenever 

either a new subband with a terra corresponding to either + or - terra 

L L 


zn ^ zn 

comes in. The results for the exact numerical integrations for — 7 — -r are shown 

T( V 

in Figure 4(b) and the curves look smoother than those in Figure 4(a). (The 
small fluctuations in the curves of Figure 4(b) are due to numerical errors in 
the multiple integrations.) The electrons in the highest subband require the 
lowest threshold energy for impact ionization and there are more incident hot 
electrons with enough energy to impact ionize these electrons. 


4.2. The Effect of the Quantum-Well Width L . 
z 

In Figures 5 (a) and 5(b), we decrease L from 200A to 80A keeping the doping 

z 

18 —3 

concentration the same (5.0 x 10 cm ). The upper bounds of the ionized 
rates are shown in Figure 5(a) and the exact numerical results are shown in 
Figure 5(b). There are three subbands for this structure and only the lowest 
one is occupied at this doping level. The ionization rate is smaller than the 
previous case in Figure 4(b). Roughly speaking, it is because the number of 
electrons per unit area in the quantum well region, , is smaller when is 

decreased and kept the same. Although the ground level of a thinner quantum 
well is higher than that of a thicker well in Figure 4(b), the higher subbands 
in Figure 4(b) also contribute to the total ionization rate and thus the ioniza- 
tion rate in the case in Figure 4(b) is larger. 


4.3. The Effect of the Band Edge Discontinuity . 

In Figures 6 (a) and 6 (b), we show the upper bounds of the ionization rates 
for all subbands and the results of the numerical integrations by decreasing the 
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A£ mole fraction x from 0.35 to 0.25 with L the same as that in Figures 4(a) 

z 

26 27 

and 4(b). We use the 60% - 40% rule * for the conduction and valence band 
edge discontinuities in our calculation. If one uses a different rule with a 
number between 60% to 85% for the conduction band edge discontinuity, the ioni- 
zation rate will decrease because the well is deeper. Once AE is determined, 

c 

the threshold energy for impact ionization can be obtained. Decreasing x will 
decrease AE^, thus the threshold energy is lowered and the impact ionization 
rate is increased. This is clearly shown in Figures 6(a) and 6(b) as compared 
to Figures 4(a) and 4(b). 

4.4. The Effect of Doping and the Total Ionization Rate . 

If one increases the doping concentration in the quantum wells, one expects 
that the ionization rate will increase because there are more electrons available 
for impace ionization. This is shown in Figure 7 where the dashed line is for 
the upper bound of the ionization rate and the solid line is from exact numeri- 
cally integration. Increasing will increase the number of occupied subbands 
also since the Fermi level will be raised higher. For each N^, the computation 
takes about 1500 cpu seconds as an average in the Cyber 175 Computer using about 

c f 

1.6 x 10° number of dicretized points (40 for variable > 40 for k^, 20 for 
4> 1 2 and 50 for k^ ) • When the Fermi level is raised, there will be more elec- 
trons in the higher subbands also and they are easier to be ionized since the 
minimum threshold energy is decreased. The other effect is that a higher doping 
concentration may also increase the tunneling probability of the bound electrons. 
This may also enhance the ionization rate. 
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5. POTENTIAL FOR DEVICE APPLICATIONS 

Avalanche photodiodes which in their operational characteristic come close 

2—fi 

to the ideal of the conventional photomultiplier have been proposed by Capasso 
and have been termed solid state photomultipliers. The principle of these 
devices is based on a highly asymmetric ionization coefficient for electrons (a) 
and holes (0). A conventional photomultiplier does, of course, not involve any 
holes at all. The question arises therefore if the new technological possibili- 
ties which are offered by MOCVD and MBE crystal growth techniques can be used 
to construct the complete solid state analogy of a photomultiplier. This ques- 
tion was one of the driving forces of our calculation in the previous sections. 
Instead of photomultiplier electrodes one would use quantum wells and the second- 
ary electrons would be emitted into the neighboring A£GaAs instead of into 

vacuum. The results for the ionization rate in Figures 4-7 show indeed that 

19 -3 

appreciable gain can be achieved if the doping of the quantum well is ~10 cm 
and the number of layers is larger than ~100. From Figure 7, the average ioniza- 

7 in 

tion rate varies from 10 to 10 l/sec. Thus the electron ionization c'oeffi- 
dents may vary from around 1 to 10 /cm or higher assuming the average electron 
velocity is 10 cm/sec or lower. The ionization rate can also be increased by 
decreasing the band edge discontinuity. Primary electrons can be generated by 
the absorption of infrared radiation either by the electrons in the valance band 
(interband absorption as shown in Figure 1) or by the quantum well electrons 
which make a transition from bound to propagating states. The latter type of 
free carrier absorption is of first order (the quantum well provides the momen- 
tum and phonons are not necessary) and therefore significantly stronger than the 
free carrier absorption in bulk semiconductors. This fact has been noted by 
Chiu et al. who also show that multi-quantum well structures (similar to our 


- 21 


structure in Figure 1, except that they also dope the GaAAAs regions) can be 

n 

used as infrared detectors . Chiu et al. did not, however, investigate the 
possibility of gain by the avalanche process. But, the possibility of varying 
the band-edge discontinuity (and therefore the range of infrared detection 
almost at will in the range of energies 0 _< hoa £ 0.4 ev) and the possibility of 
gain (analogous to the conventional photomultiplier) are enormously attractive 
features. 

There are, however, also enormous technological problems which need to be 
overcome to realize such a device. Since holes are not involved, the electrons 
which are ionized out of the well need to be replenished. This means that 
separate ohmic contacts have to be applied to a number of sets of quantum wells 
without changing the high resistivity of the AJIGaAs layers in between. It may 
be possible to accomplish -this complicated task with selective Schottky barrier 
contacts which would be ohmic for the highly doped quantum wells. Another 
severe problem is the require ment of highly doped quantum-wells neighboring to 
essentially depleted and highly resistive AJIGaAs. Both problems are not basic 
in nature and do not seem out of reach considering the recent progress in micro- 
structure technology. However, more studies (including also the effects of 
field ionization) are needed to confirm and show the device potential of such 
structures. 

Independently Capasso has recently performed experiments and has obtained 
encouraging evidence for impact ionization of the kind discussed in this paper. 
These experiments may call for an extension of our theory to include effects 
such as the tunneling-supported impact ionization * due to the band bending 
as shown in Figure 1. The bound electrons in the quantum well may acquire a 
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small amount of energy by collision with an incident hot electron and subsequently 
tunnel into the conduction band in the AJIGaAs region. This process reduces the 
ionization threshold energy in each subband. Thus the ionization rates as for 
example shown in Figure 4(b) may start from a lower threshold value for each 
subband and may also be enhanced because of the tunnel-impact ionization process. 
More work needs to be done to predict the precise performance of such a structure. 
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6. CONCLUSIONS 

The impact ionization across the conduction band-edge discontinuity of 
quantum-well heterostructures has been studied. A geometrical configuration 
to illustrate the contribution to the electron-electron interaction in the 
momentum space has been shown. The average ionization rate has also been calcu- 
lated for various doping concentrations and geometries. We have included the 
contributions from the multi-subbands and the nonparabolicity effects in our 
calculations. The application of this impact ionization mechanism to a solid 
state superlattice photomultiplier has also been discussed. 
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APPENDIX A 

THE SCREENING PARAMETER q 


The dielectric screening parameter in the static limit is given by 


[31] 


2 . -0 i /2 

i - [f- / *< E > 7 r> dE l 


(Al) 


which is the Thoraas-Ferrai screening parameter, and it reduces to l/(Debye length) 
for a Maxwellian-Boltzmann distribution in the case of a nondengenerate electron 
gas. 


When the electron concentration is degenerate, and the temperature is low, 
we obtain the Thomas-Fermi screening parameter using a three-dimensional density 
of states. 


.m*e 2 1/2 r3n"]l/6 

q - bd 


where n is the electron concentration. 


(A2) 


If we consider a two-dimensional density of state § 20 ^^’ we can carry out 
the integral in (Al) since § 20 ^^ * s a constant an< * obtain 


(^f) 

e*t 


*e 2 1/2 , N 1/2 


z 


(A3) 


where N is the number of occupied subbands. Our calculations for a quantum well 


of size L = 200 A and AE s 0.3 eV show that Eq. (A2) and (A3) give almost the 
z c 


same result for a large n. 
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APPENDIX B 


THE QUANTIZED ENERGY LEVELS E IN A QUANTUM WELL 
INCLUDING THE EFFECTIVE MASS DIFFERENCE 


The energy level E are obtained by using the boundary conditions 

♦a = *b 


^JL 1 8 * 

* 3z 


m 


a l 'b 

* 3z 

m. 


(Bl) 

(B2) 


at the two interfaces z = 0 and z ■ L z » The result is the well-known eigenequa- 
ation: 


tan kL = 
z 


2a 


1 - a 


(B3) 


where 


* 1/2 

k - (2 hlE) 7 /It 

D 


(B4) 


m* AE c - E 1/2 
° = ( i ) 


(B5) 


m 


The eigenequation (B3) includes both the even and odd solutions for the wave 
functions. The number of quantized levels in the quantum well is easily deter- 
mined from the criterion: 


If 


* 1/2 

(N-I)tt ( 2m b *E c ) l z Nw 
2 ^ 2fi ^ 2 ' 


(B6) 


the number of quantized levels » N 
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APPENDIX C 

* 

AVERAGE EFFECT MASS IN GaAs REGION (REGION b) WITH 
THE NONPARABOLICITY EFFECT INCLUDED 


The E-k relation including the nonparabolic effect is 
- 2.2 

E (1 + oE) = (Cl) 

2m 

-1 * 

where we use a = 0.67 (eV) , m = 0.063 at 77° K, and ra^ is the electron 
mass In free space. 


Using the definition 



* 2 2 
m eff (E) 9k 


and 



AE 


m e f f ( E ) dE 

~AE 


(C2) 


(C3) 


where AE^ is the conduction band edge discontinuity, we obtain analytically 
* * 2 2 

m b = m (1 + aAE c > (1 + 2aAE c + 2a AE^) (C4) 


The above procedure accounts for the nonparabolic effect up to the energy AE . 

* 

We then use this average value m^ in the parabolic band structure for GaAs in 
our calculation for the ionization rates. 
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FIGURE CAPTIONS 

Figure 1. Geometry of a structure exhibiting impact ionization across the 
band-edge discontinuity of quantum-wells. 


i 

Figure 2. The initial and final states for the electron-electron 
interaction in a quantum well structure. 


i 

Figure 3. 


A geometrical configuration in ic-space to illustrate the electron- 
electron interaction integral. With both the energy and the 
momentum conserved, the four electron wave vectors are on the same 
spherical surface. The integrations over k^ and are restricted 


to the shaded portion with the z-component between k^ and k^ + 
zn 


Figure 4(a). The upper bound of the impact ionization rates ^ as functions 

0 u 

of the incident hot electron energy E. (s E. ) for each subband and 

0 Oz 

the sum of all three subbands. The parameters used are L = 200 A, 

z 

L = 1200 A, N = 5 x 10 18 cm~ 3 , x = 0.35, T = 77°K, T = 2300°K (m* = 

D ’ ' * e a 


0.09m Q , m^ = 0.105m 0> AE c = 0.262 eV). 


Figure 4(b). The exact impact ionization rates from numerical integration of Eq. 

(25) and (42). The same material parameters as those in 4(a) are 
used here. 

Figure 5(a). The upper bound of the impact ionization rate as a function of the 

incident hot electron energy for a different quantum well width 

18 -3 * 

L = 80 A. (L = 1200A, x = 0.35, N = 5 x 10 cm m = 0.09m n , 
z D 9 a u 

* 

m = 0,105m » AE 35 0,262 eV), There is only one occupied subband 


at this doping level 
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Figure 5(b). The exact impact ionization rate as a function of the incident hot 
electron energy E^ obtained from numerical integration of Eq. (25) 
using Eq. (42). The same parameters as those in Figure 5(a) are 
used here. 


Figure 6(a). The upper bound of the impact ionization rates as functions of the 

incident electron energy E^ for a smaller Barrier height (x = 0.25, 

* * 

AE = 0.187 eV) than that in Figure 4. We use m = 0.105 m_, ra = 
c a U b 

0.083 nr for this value of x and AE . The other parameters are L 
0 c z 

1 ft — q 

200A, L * 1200A, N d = 5 x 10 cm . 


Figure 6(b). 


The impact ionization rates as functions of the incident electron 
energy E^. Same parameters are used as in Figure 6(a). 


Figure 7. 


The average ionization rate <— > (solid curve) and its upper bound 

<— > (dashed curve) as functions of the doping concentration N_. in 
t u D 

in the quantum well. We use L = 200 A, L » 1200 A, x = 0.25, T = 

z 

* * 

77° K, T = 2300° K. The parameters m , ra, and AE are the same 
e a * b c 

as in Fig. 6(a). 
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